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PREFACE 


This report is Supplement No. 1 to the Final Report (Part I) 
under Contract No. F19628-73-C-0130. This report presents the results 


of computations of cross sections and angular scattering data for 
homogeneous and inhomogeneous particles using Mie theory. This report 
gives the results of Mie calculations for single particles as well as 
macroscopic scattering data for different aerosol size distributions 


in the visible and near infrared. 


Supplement No. 2 of the Final Report (Part I) tabulates 


scattering data for different types of aerosol particles in the infra- 
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red using a r size distribuvion. 
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I. INTRODUCTION 


Light traneport in the atmosphere is remarkably influenced 
by aerosol particles. The optical effects of aerosols depend upon 
their physical structure, their chemical composition (index of re- 
fraction), their size distribution and concentration. The dependence 
of the sky intensity and its degree of polarization on turbidity 
and the size distribution of aerosol particles has been shown by 


numerous authors, Refs. 1 and 2 are excellent examples. 


The interaction of light with aerosols may be computed by 
Mie theory if one considers the aerosols as spherical particles. 
Although this may not be true for a single particle, it may be 
assumed that the aerosol particles are randomly distributed; the 
assumption of aerosol particles being spheres is generally accepted. 
Mie's theory (Ref. 3) originally assumed homogeneous particles and 
was extended later on to particles of concentric spheres (Refs. 4 and 
5) and to particles with the index of refraction varying continuously 
with the radial distance from the particle center (Ref. 6). The 
scattering properties of homogeneous aerosols and of particulates 
consisting of two concentric spheres are of utmost importance for 
the atmosphere. Currently, there exist four RRA programs that compute 
scattering data for spherical particles. These programs, designated 
as MIE-2, MIE-3, MIE-4, and MIE-5, compute cross sections and angular 
scattering data for homogeneous particles and "two-layer particles" 
consisting of a spherical kernel and a spherical shell of different 
material. All four versions of the Mie programs allow for integrating 
microscopic data over a given size distribution. The programs are 
based on the fundamental theory outlined by MIE (Ref. 3), Deirmendjian 
(Ref. 7), Aden and Kerker (Ref. 4) and Fenn and Oser (Ref. 8). A 
detailed description of the codes is given in Ref. 9. In addition 


to these four programs, a new program, designated as MIE-6, has been 


+ nee 
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written in order to compute scattering functions by combining the 
scattering functions for different size distributions and for 


different indices of refraction. 


Sample calculations with these programs are given in this 
report for inhomogeneous particles (Sections If and III), aerosol 
size distribution in the troposphere (Section IV) and stratosphere 
(Section V), and fog particles (Section VI). 
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II. EFFICIENCY FACTORS FOR HOMOGENEOUS 
AND INHOMOGENEOUS PARTICLES 


re 
| 

L 

4 The MIE-5 program was used for calculztions of the efficiency 
: factors (see Ref. 9 for definition) for particlzs of different physical 
and chemical composition. The results are preeented in Figs. 1 

through 10 as function of the size parameter, x = 2mr/\, where r is 

the particle radius and A the wavelength of the radiation incident 

on the particl:. In the graphs, the indices of refraction of the 

inner material (kernel) of the particles are designa.*¢ by MI and the 
refractive indices of the shell material by MS. 


Some of the indices of refraction cnosen for the computa- 
tions are based upon wavelength~-dependent measurements of the real 
and imaginary part of the refractive index for aerosol particles 
taken by Hanel (Ref. 10) and Fischer (Ref. 11). The index of re~- 
fractioa for “average aerosol in summer 1966 of Mainz, Germany" is 
given in Ref. 12 as 1.68 - 0.04014 for 0.3u wavelength light and as 
1.61 - 6.0721 for } = 1.0u. Figure 1 shows the well-known shape of 
the curves giving the efficiency factore of homogeneous particles 
as function of the size parameter (see Kef. 13, for example) with 
peak values of the efficiency factors for extinction occurring for 
size parameters between 1 and 10 and with a limiting value of 2 for 
large particles. In the example given, the efficiency factor for 
absorption peaks at x = 20 for m = 1.68 - 0.041 and A = 0.3 and at x = 7 
for m = 1.61 - 0.072i and dX = 1.0u where the peak values are slightly 
higher than 1. The refractive indices used for the computations given in 
Fig. 1 are for dry aerosol. In order to describe realistic conditions, 
the growth of the aerosol particles with the relative humidity and the 


related changes of the index of refraction are of special interest. 


The ratio of the particle radius at a given humidity f to ftl.at at 
f = 0 for "Mainz aerosol" is given by Hanel (Ref. 14) to be 


ee he cere yn TERT 

ae a . _— FC GS FRR BE * 
ny page ET PRS EY DL Ee PIE PG FOOLS I EE 

woe GETS ERIE TSE ETE FT EES aed a 


EXTINCT CON 


| 


ROSORPT ION 


EFFICIENCY FACTOR 
oe 


NX 
“y= 0.3u 


i 


10° 1° 107 10? 
SIZE PARANETER 


Fig. 1. Efficiency Factors for Dry Aerosols (MI = MS = 
1.680 - 0.0401 for X} = 0.3u and MI = MS = 
1.616 - 0.0721 for X = 1.0u 
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4 Fig. 2. Efficiency Factors for Inhomogeneous Aerosol Particles at a 
q Relative Humidity of 75% (MI = 1.680 - 0.0401 and MS = 
q 1.358 - 0.0001 for A = 0.3 and MI = 1.610 - 0.0721 and 

MS = 1.324 - 0.0001 for A = 1.0y) 
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Fig. 3. Efficiency Factors for Homogeneous and Inhomogeneous 
Aerosol Particles at a Relative Humidity of 75%: 
A=0.3u (MI=MS~1.489 - 0.015I for Homogeneous Aerosols 
and MI=1.680 ~- G.0401 and MS=1.358 - 0.COOI for 
Inhomogeneous Aerosols) 
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Fig. 4. Efficiency Factors for Homogeneous and Inhomogenesus 
Aerosol Particles at a Relative Humidity of 75%: 
Aw1.0U (MI=MS*1.441 - 0.0291 for Homogeneous Aerosols 
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et ioe OR IP IIIT RETIN 9 AEs TET ERE EST Te 
GET LTTE TCP TTITE ET FET PR =F or ae 


a ae ve TT FE EE EN 


6 


8 anaes Particles 


ABSOxe TIO 


Homogeneous Particles 


EFFICIENCY FRCTOR 
° 


Inhomogeneous wie 
—~RH=99.9% 

7 4 

1c! 1a! 10 10° 
Sftzt PARAMETER 


Fig. 5. Efficiency Factors for Aerosols at 95% and 99.9% 
Relative Humidity: A=1.0y (When RH=95%, MI=1.620 - 
0.0721 and MS¥1.324-0.000I for Inhomogeneous 
Particles and MI=MS=1.360 - 0.009I for Homogeneous 
Particles. When RH=99.9%, MI=1.620 ~ 0.0721 and 
MS=1.324 - 6.0001 for Trhomogeneous Particles and 
MI=MS=1.325 - 0.000I for Homogeneous Particles) 
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Fig. 7. Efficiency Factors for Water Particles Coated 
With Material of Higher Refractive Index: 
r/r, = 1.348 
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Water-Suspenced Carbon Particles 
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Fig. 9. Efficiency Factors for Hydrosol Particles 
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Fig. 10. Eificiency Factors for Water Bubbles 


r(f)/r, = 1.348, 2.0 and ~ 8.0 for relative humidities of £ = 75%, 
95%, and 99.9%, respe-tively (the ratio r(£)/r, ks somewhat dependent 
upon the particle size for high humidities which will not be con- 
sidered in this study). The average index of refraction at humidity 


f is then given by 


r(f) 
a ea Oo ae, ae. 


where a, and mn, denote the indices of refraction for water and dry 
aerosol. The efficiency factors for "Mainz aerosol" are given in 
Figs. 2, 3 and 4 for £ = 75% and in Fig. 5 for relative humidities 

of 95% and 99.9%. For the data shown in Figs. 2, 3, 4 and 5 for 
inhomogeneous particles, the index of refraction of the kernel. was 
taken as that of the dry aerosol; the shell material was assumed 

to be water with the ratio of the particle size to the size of the 
kernel chosen as r(£)/r). The data shown in Figs. 2, 3, 4 and 5 

for homogeneous particles were computed using the average index 

of refraction given by Eq. (1). For f£ = 75%, differences 

in the calculations are seen mainly for size parameters 

smaller than ~ 30 with the efficiency factors for the homogeneous 
particies being smaller for absorption and J :rger for extinction than 
those calculated for the spherical-shell particles. For f = 95% 

(Fig. 5), some differences in the dependency of the efficiency factor 
for extinction upon the size parameter are seen for all size para-~ 
meters x > 3. Homogeneous parcicles are seen in Fig. 5 to be more 
absorbing chan shell particles, especially for large particles with 
the absorption efficiency factors for x = 125.66 calculated to be 0.95436 
and 0.43785 for homogeneous and inhomogeneous particles, respectively. 


For f = 99.9% (Fig. 5), no apparent differences were found between the 


efficiency factor calculations for both homogeneous and inhomogeneous 


particles in the ent‘re range 2.5133 < x < 113.0973. 
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The efficiency factors for homogeneous carbon particles 
and for carbon particles surrounded by water shelis with the ratio 
of the maximum and minimum shell radii taken to 1.35 are shown in 
Fig. 6. The graph shows the well-known smooth curves without the 
oscillations for the pure carbon particles (see also Ref. 8). The 
cu ves shown in Fig. 6 for the shell particles make it apparent 
that the observed oscillations in the extinction efficiency factor 


are caused by scattering and not by absorption. 


Figure 7 gives the efficiency factore for water particles 
coated with a material of higher refractive index, witi the ratio 
of the particle to the kernel size taken to 1.348. These data 
were previously reported in Ref. 15, pp. 90/91, with the indices 


of refraction erroneously given in reversed order. 


For particles suspended in media other than gases, the 
effective index of refrartion is obtained as the ratio of the 
refractive indices of the particle, a,» and of the surrounding 


medium, m_: 
n 


nm, =z a) = (, -_ ik) /(a, = ik) = 


aa +kk -i@k -ak) 
m m m Tm 


. (2) 


For the case of carbon particles in water and 0.3 wavelength light, 
Eq. (2) leads to 


1.95 0.66 


e * 1.358 ~ 1.359 4 = 1-436 - 0.4864 : 
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Examples of the scattering and extinction by such hydrosols are 
given in Figs. 8 and 9. Figure 8 sefers tc gas bubbles in water 
(the refractive index of 0.775 corresponds to a waveiength of 

approximately 2.25} and to carbon particles; Fig. 9 to “aerosol 
particles" in water. ‘The indices of refraction shown in Figs. § 


and 9 are “effective” indices of refraction given by Eq. (2). 
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Figure 10 shows the application of the MIE-5 program to 
“water bubbles", i.e. water spheres having gaseous kernel. The 
ratio of the outer radius to the inner radius was taken to be 
1.35 for A = 0.3u and 1.01 for A = 1.0n. 
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TIZ. SCATTERING PROPERYLES OF HUMID AEROSOL PARTICLES 


The data in Section II revealed differences in the efficiency 
factors for homogeneous and inhomogeneous particles. In particular, 
it was found that there are differences in the efficiency factors 


for aerosol particles consisting of a spherical kernel of water- 


insoluble material and a spherical water shell and for homogeneous 
particles assuming an average index of refraction that is calculated 
by mixing rules, It is of specia] interest how macroscopic scatterimp 


data would compare when calculated with these two methods. 


The data given in Figs. 2 and 3, corresponding to a relative 
humidity of 75%, were integrated over a power-law size distribution, 


n(c) = 1.0628r74 3 


with the limiting radii taken as Tada 0.05392u and Pr ee 13.48u. 
For the spherical-shell model the ratio of the outer and inner radii 
was assumed to be 1.348, The ratio of the mass of the aeros 9] 
particle at x % relative humidity to the mass at 0.0% relative 


humidity is ziven by 


In p, 3 
W 
— ==] 4— xe ~ 1 
m po r 
ke) ° 


P 
7 


; 2 
where e and a, are che densities of water (0. = 1 g/cm") and the dry 
aerosol particles, respectively and rir ig the ratio of the outer to 
the inner radii of the spherical shell particles. For t/r, = 1.348 


and a relative humidity of 75%, the mass ratio is given by 


i, 1.449 
Sh. yg Lea 
M4, p 
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The density of dry aerosol particles, Oo? is assumed to be 3.0 g/cm 


in Ref. 12 (Model 1). For the homogeneous model, Eq. (1) was used 


to calculate the index of refraction. The computations were made for 


5 wavelengths ranging from 0.3 to 1.0u, the indices of refraction 


for dry aerosols were taken from Ref. 12 (model 1). 


The calculated optical parameters from these calculations 


are listed in Table I. The tables lists the indices of refraction, 


the actenuation coefficients (pyr ¢ denote the ex- 


> Oscap? and yp 
tinction, scattering, and absorption coefficients, respectively) , 


the single scattering albedo, the average cosine of the scattering 


angle, cos0, and the maximum value of the degree of polarization, 


Pax’ As expected from the results presented in Figs. 2 and 3, the 


absorption coefficients are higher for the spherical-shell particles 


when compared with those for the homogeneous particles with an 


average index of refraction. The scattering coefficients are found 


to be smaller for the spherical-shell particles. This produced 


somewhat smaller extinction coefficients and smalJer values of 


the single scattecing > edo. With the exception of 0.3) wave~ 
length light, che maximm degree of polarization is smaller for the 


inhomogeneos particles for all wavelengths being considered. 


The normalized phase matri.ces are given for all wavelengths 


in Tables Al through A5 for spherical--shell particles and in 


Tables A6 through ALO for homogeneous particles. The scattering and 


prlarizatio:r functions for 0.3y and 1.0 wavelength light are pre- 


sented in graphical form in Ref. 15 for a similar problem. It was 


found in Ref. 15 that larger differences in both computations 


occur f¢.: both the scattering function and the degree of polariza~ 


tion, especially for scattering angles larger than 100 degrees. Ths 


effects of inhomogenuity of the shell particles do not appear to be 


important for scattering angles less than about 4°. The scattering 


function for the homogeneous particles exceeds those for the shell 


particles for scattering angles between 4° and 40° by up to 20%. At 


scattering angles greater than about 45° the scattering function for 


shell particles exceede that computed for homogeneous particles. 
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Figure 11 presents plots of the ellipticity, tanf, as function 
of the scattering angle. For the data in Fig. J1 it was assumed that 
the radiation incident on the particles is linearly polarized with 
the plane of polarization inclined to the reference plane by 45°. 

The ellipticity, which is dependent on the four elements of the phase 
matrix, is given by the equation 


~i 


tang = ——-—- : (3) 
/ 2 2 2 / 2 2 
0.25(4,-1,) + i, + i, + 0.25(i,-1,) + i, 


It was found that the ellipticity is more sensitive to the homo- 
genuity of the particles than are the scattering function and the 
degree of polarization (Ref. 15). The data in Fig. 11 show a shift 
of the peak value of the ellipticity for 0.3u wavelength light from 


122° to 160° scattering angle when two-layer particles are assumed. 


The data presented in Sections II and III show large differences 
in the scattering properties, especially in the yolarization features, 
for “Mainz aerosol" at a relative humidity of 75% when homogeneous and 


inhomogeneous particles are assumed. The differences are expected to be 


even larger for a relative humidity of 95% (see Fig. 5). Since natural 
aerosols consist of roughly 50% water insoluble material (Refs. 16 and 
17, the ratio of water~solubles to water insoluble materials shows 
large fluctuations and is highly dependent upon the origin of the 
particles), care has to be taken in interpreting data of the elliptical 
polarization (Ref. 18). Previous studies of the scattering properties 
of aerosol particles as function of the relative humidity, such as those 
reported in Ref. 12, should be extended to consider the scattering 


properties of water-insoluble materials as well as water solubles. 
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IV. SCATTERING PROPERTIES OF TROPOSPHERIC 
AEROSOL PARTICLET 


The size distribution of continental aerosol particles in the 
troposphere may be classified into three different scattering types, 
designated in Ref. 19 as types A, B, and C. According to Ref. 19, 
different scattering types result mainly from different size distributions 
of the aerosol particles in the region 0.15u < r < 0.60, where r is the 
particle radius. Size distributions of scattering types A, B, and C 
are obtained through a concave, linear, and convex shape of the particle 
size distribution (on double logarithmic scale), resulting generally ir 
a concave, linear and convex shape of the curve giving tke relation 
between extinction coefficient and wavelength, and also resulting in 
increasing, constant, and decreasing values of the scattered intensity 
ratios, I(red)/I(blue), for decreasing scattering angles in the range 
from 50° to 10°. Data are given in this section for the angular 
scattering of aerosol particles having different size distributions. 
Though listed under a specific scattering type, some of the size dis- 
tributions used are transitions between different scattering types 
(type AB, etc.). 


4.1 Scattering Functions for Type B Aerosols 


Type B aerosols are described by power-law size distributions. 
Scattering functions for aerosols with power-law size distributions 
are published by numerous authors, see for example Ref. 20. Most of 
these data, however, are obtained assuming that the absorptiou of light 
by aerosol particles may be neglected, an assumption which, according to 


earlier mentioned measurements by Fischer (Ref. 11), is very doubtful. 


Some model calculations were therefore accomplished for type B aerosol’ 
and indices of refraction of 1.50 - 0.01i and 1.55 - 0.01i (similar 
refractive indices are obtained for models 5 and 6 in Ref. 12 for relative 
humidities of f ~ 50%). 


t 
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For the calculations, the parameter v in the size distribu- 
tion 


eee e (4) 


was taken as 4.0, 4.5, and 5.0. The parameter c was taken such that 
the particle number was constant for each size distribution at a 
particle radius of 0.55 micron. The lower limit of the particle 
radii was varied between 0.0ly and 0.12, the upper limit was taken 
to 10u or 20H. The tctal number of particles per unit volume is 


given by 
- r 
se ag) eee 2 (-v+t1) (-v4+1) 
N=c } Per | aes aaa Tee 7 (5) 
x 
min 


A complete description of the models used is given in Table IT. 

Values of the attenuation coefficients, as obtained from the MIE~2 
program, are listed in Table III; the corresponding scattering functions 
are tabulated in Tables Bl through B28. 


Tne scattering functions listed in Tables Bl through B28 were 
used to calculate the sky radiance in the solar almucantar for each 
of the models being considered; the computations were performed with 
the FLASH and BRITE Monte Carlo Programs (Ref. 15). The aerosol 
optical thickness was chosen such that the optical thickness would 


equal 0.23 for the size distribution with v = 4.5, fai 0.06um and 


for 0.55u wavelength light. The aerosol optical ae for other 
models and wavelengths is then given by multiplying 0.23 with the 
ratio of the extinction coefficient for that model and the one for 
0.55u, v= 4.5 and Fes = 0.06um. The atmosphere was assumed to he 
homogeneous, the angle of incidence was taken to be 70°. The receiver 


azimuth angle, a, is then related to the scattering angle, $, through 


: { 
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LABLE III. ATTENUATION DATA FOR TYPE B AKROSOLS 
m= 1.50 - 0.01i 


Ne ne ee a ee eed 


-1 
m v Tin) ACh) Sooay (cm ) Onyg Com ) 
1.50-0.01i 4.0 0.06 0.30 7.81-10 8.58-10 
0.45 5.52--10 6.07-10 
0.70 3.61-10 3.98-10 
1.60 1.57-10 1.74-10 
1.59-0.014 4.5 0.06 0.30 1.45~09 1.56-09 
0.45 8.75-10 9.44-10 
0.70 4.69-10 5.12-10 
1.60 1.38-10 1.54~10 
1.50-0.01i 5.0 0.06 0.30 2.94-09 3.13-09 
0.45 1.54~09 1.65-09 
0.70 6.87-09 7.52-10 
1.60 1.37-10 1.60-10 
1.50-9. 014 4.5 0.01 0.30 1.71-09 1.92~09 
0. 45 9 e 29-10 1 e 06-09 
0.70 4.79-10 5.61-10 
1.50-0.01i 4.5 0.12 0.3¢ 8.48-10 9.20-10 
0.45 6.83-10 7.34-10 
0.70 4.27-10 4.59-10 
1.60 1.36-10 1.49-10 
1.55-0.014 4.5 0.06 0.45 9.75-10 1.05-09 
0.55 741-10 7.98-10 
1.60 1.55-10 1.71-10 
1.55-0.0114 5.0 0.06 0.45 1.75-09 1.87-09 
0.55 1.23-09 1.32-09 
0.70 7.86-10 8.53-10 
1.60 1.57-10 1.80-10 
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cosd? = 0.1170 + 0.8830 * cosa . (6) 


Some of the Monte Carlo calculated data have been reported 
in Ref. 15, Out of all of the models considered, it was found that 
the model with v = 5.0, a 0.06ym and an index of refraction of 
1.55 - 0.01i (model 504) gives the best approximation to observed 
scattering functions and color ratios for type B aerosol (Ref. 19). 
The sky scattering functions for this model are plotted in Figs. 12 
and 13 for single scattering (SS) and for total scattering assuming 
ground albedo values of A = 0.0, 0.2, and 1.0; the calculated aerosol 
scattering function (dotted line; normalized to give the same value 
as the sky scattering function for 0.1 degrees) is given for compari- 
son. As one would expect, the sky scattering function is more affected 
by molecular and multiple scattering for 0.45y wavelength light 
(Fig. 12) than it is for 1.60u. The color ratios, 1(0.7u)/1(0.45y) 
and 1(1.6)/1(0.7u), are given in Fig. 14. 


Measurements in Europe under visual conditions usually produced 
by Type B aerosols (Ref. 21) gave results in the near IR that were not 
in agreement with the Mie calculations described in this report for 
aerosols with Type B size distributions. More recent measurements in 
Nova Scotia and Bedford, Massachusetts (Ref. 22) conform to the results 
presented here for } = 1.60um. A very low refractive index or a large 
influence of a possibly existing small deviation of the size distribution 
at r 2 lpm from the power law had to be assumed to explain the steep IR 
scattering functions in the European data. More nocel calculations are 
required to resolve the problem. 
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Fig. 14. Color Ratios in the Solar Almucantar for - Size 
Distribution (1{0.7y)/1(0.45u), upper graph; 
1(1.60u)/1(0.7u), lower graph) 
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4.2 Scattering Functions for Type A Aerosols 


Three models, denoted as models 200, 500, and 600, have been 


chosen as examples of the size distributions for aerosol particles of 


5 
; 
3 
a 
2 
$ 
4 
; 
4 
5 
é 
4 
Fy 


scattering type A. The size distributions used are tabulated in Table 
IV and plotted in Fig. 15. For the computations, a total of 600 radii 


enti Demande tere 


was used for the integration over the particle size distribution. Two 
methods were used to interpolate the number density of a particle of 
radius r within the interval r, and 7) (where r 


1 1 
given in Table IV). For model 200, the interpolation was done using 


and ty are the radii 


the formula 


log(r/r,) 


n° 2 (7) 
7 


n(r) = n(r,) + (n(x,) - nro) Tote, /r. 


The number densities for models 500 and 600 were interpo- 
lated through 


[ \* 


n(r) = n(r,) es r, ’ 
where 
log(n(r,)/n(r,)) 
Oa (8) 
log(r,/r,) 


Tae index of refraction for type A aerosols was taken to be 1.55 - 
0.01i. Attenuation coefficients for type A aerosol particles as a 
function of wavelength are given in Table V, the scattering functions 
ure listed in Tables Cl through C12. The scattering functions are 
pletted in Fig. 16 for all three models and for 0.7, wavelength 
light, the color ratios (1(0.7u)/1(0.45y)) are given in Fig. 17 


(the color ratios were obtained from the normalized phase-functions 
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TABLE IV. SIZE DISTRIBUTIONS FOR TYPE A AEROSOL MODELS 


ul Model 200 Model 500 Model 600 
.06 99.9 102.6 940.3 
.07 72:2 74.1 582.5 
.08 43.6 44.8 307.8 
.09 30.4 31.2 190.6 
ee 22.3 22.9 126.1 
15 2.88 2.95 10.8 
12 .640 0.657 1.81 
.25 .223 0.229 6.3505 
3 091 0.090 0.171 
4 .025 0.028 3.53-2 
5 O11 0.012 1.23-2 
.6 6.20-3 5.69-3 584-3 
4 3.93-3 2.75-3 3.17-3 
8 2.70-3 1.60-3 1.91-3 
.9 1.97-3 1.01-3 1.23-3 

i 1.42-3 7.164 8.04-4 

1.5 5.13-3 2.04-4 1.93-4 

2. 2.77-3 8.99-5 7.82-5 

2.5 1.79-4 483-5 4.04-5 
3. 1.31-4 3.21-5 2.47-5 
4.0 7,79-5 1.43-5 1.10-5 
5. 5.09~5 7.71-6 5.75-6 
6. 3.59-5 483-6 3.38-6 
Vs 2.66-5 3.21-6 2.15-6 
8. 2.15-5 2.29-5 1.51-6 
9. 1.79-5 1.84-6 1.12-6 

1.43-5 1.40-6 8.07-7 

7.22-6 400-7 2.72-7 


4.92-6 1.50-7 1,39-7 
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Fig. 15. Model Size Distributions for Type A Aerosols 
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Fig. 16. Normalized Scattering Functions for Type A Aerosols, 
A = 0.7uUpy 
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Fig. 17. Color Ratios for Type A Aerosol Size Distributions 
(1(0.7u)/1(0.45p)) 
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and not from the volume scattering functions). It is seen that the 
forward scattering is more pronounced for the models that have 
relatively more large particles. The scattering function 1s therefore 
highest for model 200 and lowest for model 600, an opposite effect 

is observed for the color ratio (1(9.71)/1(0.45u)) at small scattering 


angles. The color ratios show peak values at scattering angles of 


LE TTR EL OF ET AATF a ARE 


~ 1° and 180°; the ratios are low for scattering angles of 0° and 


~ 40° for all models being considered. 


4.3 Scattering Functions for Type C Aerosols 


The index of refraction for type C aerosols was taken to 
be 1.55 - 0.01i. Size distributions for type C aerosol particles 
are given in Table VI and in Fig. 18 (Eq. 7 was used to interpolate 
the size distributions for models 100 and 300, resulting in curves 
that are not as smooth as those shown in Fig. 15). Attenuation data 
are listed in Table VII, the corresponding scattering functions in 
Tables D1 through D20. 


The normalized scattering functions for 0.70u wavelength 
light are plotted in Fig. 19, and the color ratios (1(0.71)/1(0.45p)) 
are given in Fig. 20. Figure 19 shows the high dependency of the 
magnitude of the forward scattering upon the (relative) number of 
large particles. All size distributions show a relative maximum of 
the scattering function at about 20° scattering angle. The scattering 
functions are almost identical for larger scattering angles. Figure 20 
shows that maxima of the color ratios (1(9.7y)/1(0.45u)) occur at 
scattering angles of ~ 1° and ~ 50° - 60°, the second maximum corres- 
ponding to the so-called Bishop-ring that is observed for type C 
aerosol size distributions (Ref. 19). 

Model 100 was used for calculating the sky radiances in the 


solar alumcantar (see also Section 4.1 and Ref. 15). Piots of the 


sky scattering functions are given in Figs. 21, 22, and 23 for 
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TABLE VI. SIZE DISTRIBUTIONS FOR TYPE C AEROSOL MODELS 


tm] Model 100 Model 300 Model 400 Model 700 Model 800 
06 9.73 9.73 30.0 89.2 89.2 
07 8.07 8.97 26.0 63.4 63.4 
-08 6.79 6.79 21.0 46.7 46.7 
-09 5.69 5.69 18.0 34.8 34.8 
ol 4.83 4.83 16.0 26.6 26.6 
015 2.57 2.57 9.0 9.42 9.42 
02 1.42 1.42 6.0 3.91 3.91 
225 -716 716 3.5 1.57 1.57 
3 - 363 - 363 1.0 0.666 0.666 
04 4.86-2 4.86-2 4.90-2 6.66-2 6.66-2 
05 1.34-2 1.34-2 1.34-2 1.47-2 1.47-2 
6 5.69-3 5.69-3 5.69-3 5.22-3 5.22-3 
7 2.75-3 2.75=3 2.75-3 2.16-3 2.16-3 
8 1.60-3 1.60-3 1.40-3 1.10-3 1.10-3 
9 1.01-3 1.01-3 8.00-4 6.17-4 6.17-4 
1. 7.16-4 5.00-4 4.50-4 3.94-4 3.94-4 
1.5 2.04--4 1.50-4 7.50-5 7.48-5 7.48-5 
2. 8.99-5 4.00-5 1.95-5 2.47-5 2.47-5 
2.5 4.83-5 1.50-5 9.00-6 1.06-5 1.06-5 
3. 3.21-5 7.00-6 4.10-6 5.89-6 5.30-6 
4. 1.43-5 2.40-6 1.40-6 1.96-6 1.60-6 
5. 7.71-6 1.10-6 5.00-7 8.48-7 7.00-7 
6. 4.83-6 5.50-7 2.50-7 4.43-7 3.40-7 
7. 3.21-6 2.10-7 1.30-7 2.57-7 1.80-7 
9. 2.29~6 1.90-7 6.06~8 1.57-7 1.00-7 
9. 1.84-6 1.30-7 3.70-8 1.12-7 6.00-8 
10. 1.42~-6 1.00-7 2.10-8 8.09-8 3.60-8 
15. 4.83-7 1.00-8 3.00-9 1.77-8 4.00-9 
20. 2.57--7 2.00-9 7.80-10 7307-9 9.00-10 
25. 2.80-10 2.80~-10 
30. 1.00-10 1.00-10 
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Fig. 18. Model Size Distributions for Type C Aerosols 
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A = 0.70p 


Eee Me ee 


eg SN SS 


ee Bey iy ee et ees 


ge Soe ee OS eR 


Rl ii Het, 


oe a id pales 


id 


10 


INTENSITY (PHOTONS/STER}) 
2, 


WAVELENGTH = 1.6000 MICRON 


ig. 23. 


ee a a ey ee ee Oe Se ey ee ey ee 


39 


RRA BRITE 
ATMOSPH. 100 


AEROSOL 


SCATTERING 
S8@® A=0.004 A=0.20 + FUNCTION -- 


10° 10! 10 
SCAT. ANGLE (DEGREES) 


Sky Scattering Function for ModeJ, 100, A = 1.60y 


ce ee Os, eee a ee, ae ee ee a ee See ede aes eee SB i oe 


Pe ae ger eee ae 


ay eS ae oe oe eee 

= ss: BN PATEL EE hah ~ hehe the eS ie. 
Deiat tee wi ag FOIE SP OE ITED TEE oO SAT IEEE OT TN eI TT oe e 
4° PERE API PRAT eS = 


we wert Sh en 


40 


A = 0.45, 0.70u, and 1.60, resnectively. The corresponding color 
ratios are given in Fig. 24. Data giving the scattering functions 

and color ratios as obtained from the MIE data for model 100 are 

shown in the graphs for comparison (broken lines, normalized tu the 
corresponding single scattering values at 0.1° scattering angle). 

It is seen that the scattering functions are more affected by molecular 
and multiple scactering at shcrter wavelengths for scattering angles 
larger than ~ 10° resulting ia a noticeable decrease of the color 

ratio in the area of the Bistop-ring. The Bishop-ring almost disappears 
for albedo values A> 0.2. ‘it is seen in Fig. 24 that the color 

ratios for Model 100 are not significantly dependent on the magnitude 
of the ground alhedo. 


A comparison with measured sky scattering functions and 
color ratios in Ref. 19 reveals that a good approximation of the 
size distribution for type C aerosols would be obtained by combining 
models 400 and 700, using model 700 for particle radii r > lum and 


increasing the number of particles for model 400 and r < ~ 0.4m. 
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Vv. SCAT ERING FUNCTIONS FOR STRATOSPHERIC AEROSOLS 


Little is known about aerosol in the stratosphere and there 
exist only a few direct measurements of aerosol size distributions 
in the stratosphere. A more recent paper by Friend (Ref. 24) shows 
that -~ at least in the 40 km layer -~ ~ 93% of the stratospheric 
aerosols are Aitken particles and - 6% are particles having radii in 
the region arcund 0.3m. Friend's size distribution for the dust 
particles can be approximated by the gamma distribution (Ref. 25) 


as 
aN, 8.5 -28.33r 
dlogr re . (9) 
or 
dN = os 7 ry 5 -28 e 33r 
= aer e (9a) 


with the peak of the distribution (9) obtained at t = 8.5/28.33 = 
0.3um and of the distribution (9a) at 0.265ym. 


Mie calculations were performed for Aitken particles 
(model 30) and for the size distribution given by Eq. (9a) (model 20). 
For the Aitken perticles, the size distribution was chosen to be 
constant (18.68) with limiting radti of r= 0.Clum and 0.06um. The 
total number of Aitken particles is therefore 18.68 * 0.05 = 0.934cem >, 
The value a in the distribution (9a) was taken to 1.0396 ° 10’ , the 
lower and upper limit of the particle radii was chosen to 0.lum and 
19.0um. The total number of particles is then approximately (Ref. 22) 


N = 1.0396 + 107 « 28,3378*> . 78.5) = 0.066cem™>. The index of re- 


fraction was taken to 1.50 - 0.0i for both model 20 and model 30. 
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The computed scattering coefficients are listed in Table VIII, 
the normalized phase matrices are given for wavelengths of 0.4in, 0.50n, 
0.55, 0.70, and 0.85 in Tables El through E5 for model 20 and in 
Tables E6 through E10 for model 30. 


TABLE VIII. SCATTERING COEFFICIENTS FOR STRATOSPHERIC AEROSOLS (cm!) 


m= 1.50 - 9.0i 


ne SS 


Wavelengths Model 20 Model 30 
0.41y 6.500-10 2.943-12 
0.50y 7.184-10 1.357-12 
0.55u 7.315-10 9.2699-13 
0.70u 6.834-10 3.548-13 


0.85y 5.747-10 1.629-13 


It can be seen from Table VIII that the scattering coefficients 
of particles in the 0.3um region is almost independent of the wavelength, 
\ whereas the scattering coefficient for the Aitken particles varies 


proportionally to iy 


(the scattering coefficient for molecules varies 
by 1"). The scattering coefficients for particles having the size 
distribution (9a) are much larger than those for the Aitken particles 
even though the Aitken particles outnumber the larger particles by a , 
factor of .934/.066 = 14.15. These high values can be verified by a 
simple hand calculation. The ratio of the scattering coefficients for 


models 20 and 30 can be approximated by 
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(oj N.. ° 0 2637 * T * Q(x,) 
220 ass, 205 2 
0 e 2 e e 
30 Nyy * 0.035" * 1 + Q(x,) 
Q(x,) 
6.05% 2 ; (10) 
Q(x,) 


where Q(x) is the efficiency factor for particles of size x = 2mr/\ 
(xr = radius, \} = wavelength). For d = 0.41) we obtain (MIE-3 program) 
Q(2m * 0.265/.41) = Q(4.06) ~ 4.05 and Q(2T7 * 0.035/.41) = Q(.536) = .02 


and therefore Foo/0 = 820; for A = 0.85y the ratio 559/% is computed 


30 
to ~ 5.03 * 103. Both of these ratios are in the same order of mag- 


nitude as those obtained from Table VIII. 


A combined phase function for modeis 20 and 30 may be calcu- 
lated by 


POO) = (pq * Pog() + O59 * Pyg($)/ (yy + F545), (11) 


where Po9 and P30 denote the phase functions for models 20 and 30, 
respectively. Ever though the phase functions for these two models 
are quite different (Fig. 25), it can be seen from Table VIII and 
Eq. (1) that the stratospheric aerosol Scattering is dominated by 
the large particles and -~ under the assumption of the particle 


concentrations for Aitken particles and large particles as reported 
in Ref. 24 -- that light scattering by Aitken particles may be 
neglected in the stratosphere. 
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VI. FOG SCATTERING IN THE VISIBLE 
AND NEAR IR 


The gamma distribution (Ref. 22) 


an = Deree (9b) 


was used for MIE calculations in fog- Four different types of fog 
particles were assumed. These size distributions are designated as 
models 1, 2, 3, and 4 (see Refs. 26 and 15 ; the constants needed 

in Eq. (9b) are listed in Table IX. Plots of the size distributions 
for fog models 1, 2, 3, and 4 are given in Fig. 24 of Ref. 15. 

Using Eq. (9b), the mean particle radius is given by 


ee (m+1) (m+2) 
2 b (12) 


and the total number of particles by (Ref. 25) 


Ny = D+ T(m+l) - pty (13) 
The extinction coefficient may then be approximated by 
o = Mme Nyt ry 10 (kn), (14) 


where No and tT, are expressed in cm and ym, respectively. 


The calculations were performed for wavelengths of 0.40n, 
0.6943 (Ruby laser wavelength), and 0.8011. The indices of refraction 
for these three wavelengths were taken to be 1.343 - 0.01, 1.330 - 
0.24 » 10° 
tion coefficients are compared in Table X with the extinction 
coefficient given by Eq. (14). The corresponding phase matrices for 


i, and 1.328 - 0.01, respectively. The computed attenua- 


all models and wavelengths under consideration are given in Tables 
Fl through F12. Plots for ) = 0.6943u and fog uodels 1, 2, 3, and 4 
are given in Fig. 26. 
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APPENDIX A 


Normalized Phase Matrices for Aerosol Particles in an 


Atmosphere with a Relative Humidity of 75% 
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APPENDIX B. 


Normalized Phase Matrices for Type B Aerosols 
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APPENDIX C. : 


= Normalized Phase Matrices for Type A Aerosols 
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APPENDIX D. 


Normalized Phase Matrices for Type C Aerosols 
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4 APPENDIX E. 


Normalized Phase Matrices for Stratospheric Aerosols 
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Normalized Phase Matrices for Fog Particles 
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